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The passage-time distribution for a spread-out quantum particle to traverse a specific region is 
calculated using a detailed quantum model for the detector involved. That model, developed and in- 
vestigated in earlier works, is based on the detected particle's enhancement of the coupling between 
a collection of spins (in a metastable state) and their environment. We treat the continuum limit of 
the model, under the assumption of the Markov property, and calculate the particle state immedi- 
ately after the first detection. An explicit example with 15 boson modes shows excellent agreement 
between the discrete model and the continuum limit. Analytical expressions for the passage-time 
distribution as well as numerical examples are presented. The precision of the measurement scheme 
is estimated and its optimization discussed. For slow particles, the precision goes like E~ 3 ^ 4 , which 
improves previous E~ l estimates, obtained with a quantum clock model. 

PACS numbers: 03.65.Xp, 03.65.Ta, 03.65.Yz, 78.20.Bh 
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I. INTRODUCTION 

Time-of-flight measurements are a standard tool for 
many experimentalists. Since the particles or atoms in- 
volved are usually fast, their center-of-mass motion is 
typically treated classically, yielding a simple description 
of the time-of-flight measurement. But as the diffraction 
and interference experiments using temporal (instead of 
spatial) slits of Szriftgiser et. al have shown such a de- 
scription of the center-of-mass motion by means of classi- 
cal physics is not always sufficient: the advance of cooling 
techniques has made it possible to create ultracold gases 
in a trap and produce very slow atoms, e.g., by opening 
the trap. Whenever such ultracold atoms are involved, 
the spatial extent and the spreading of the wave function 
can show noticeable effects. Even the seemingly simple 
question of the time spent by a particle in a given region 
of space does not possess a simple and definite answer. 
Related to this "dwell-time" problem are the problems 
of "passage time," concentrating on those particles that 
actually cross the region of interest and are not reflected, 
and "tunneling time" in the case of a barrier that clas- 
sically cannot be traversed inside the region of interest. 
These problems have on one hand been treated axiomati- 
cally aiming at ideal quantities relying only on the system 
of interest, see, e.g., Refs. @, 0, Other approaches 

may be called "operational" in the sense that a sort of 
"measurement device" is introduced to which the system 
of interest is coupled; see, e.g., Refs. @, @, @, H, OH EU • 
For a critical review on different approaches to tunneling 
times see, e.g., Refs. fl2l. [l3j. 

A distinction that can be drawn between the various 
time-related quantities recalled above concerns whether 
they are pre- or post- decoherence. Since the present 
work models the measurement apparatus, it can be 
thought of as post-decoherence; in fact it is one of 



our objectives to include decoherence- inducing processes, 
thereby eliminating some of the black magic of quantum 
measurement. On the other hand, one can take, for ex- 
ample, a path integral approach to tunneling time [l4l.[l5j 
in which the paths are sorted by (variously defined) times 
spent in the barrier region. For these paths, amplitudes 
are retained, so it is very much a pre-decoherence calcu- 
lation. As a result interpretive issues arise (discussed for 
example by [Hj]) which do not enter in the present work. 

In the present paper we investigate a particular mea- 
surement scheme for passage times, mimicking the way 
a time-of-flight experiment is typically performed: Em- 
ploying a particular spin-boson detector model for mea- 
suring quantum arrival times investigated in Ref. [l7j |. 
our measurement scheme involves two measurements of 
arrival time by means of this detector, one upon en- 
tering the region of interest and one when exiting; see 
Fig. Q] This scheme can be expected to distort the par- 
ticle's wave function less than a scheme based on the 
(semi-)continuous coupling of the particle to a measure- 
ment device. Indeed, it will be seen that a passage-time 
measurement by means of "slow" detectors yields a rather 
broad passage-time distribution. But detectors respond- 
ing quickly to the presence of the particle also yield a 
broad distribution. This will be shown to be a quan- 
tum effect involving the Heisenberg uncertainty relation. 
However, there is an intermediate range for the detec- 
tor parameters yielding optimal results. In the best case 
the precision of the measurement can be estimated to 
behave like E~ 3 ^ 4 , improving the results obtained from 
clock models QUE!. 

In Sec. [n] we review the detector model and its ap- 
plication to quantum arrival-times. In Sec. IIIII the par- 
ticle's wave function immediately after the detection is 
investigated, and in Sec. IIVI analytical formulas for the 
application of the detector model to passage times are de- 
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is the free Hamiltonian of the particle; 
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FIG. 1: The first detector detects the particle at entry into the 
region of interest, the second one at exit. From the correla- 
tion of these two arrival-time distributions, the passage-time 
distribution is obtained. The ^-positions correspond to the 
numerical examples discussed in Section [IV Bl 
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with ||) • denoting the excited state of the j th spin and 

<^IT),;(TI - ID^-UI, (4) 

is the free Hamiltonian of the detector; fuv^ is the energy 



difference between ground state and excited state of the 
j th spin, and fiwj *' > is the coupling energy between 
the spins j and k. Further, 

-Hbath = E hw t a\ai, (5) 



rived. Numerical examples are investigated in Sec. IIV Bl 
For simplicity, we concentrate there on passage times for 
free particles; the extension to passage times in the pres- 
ence of barriers, and thus to tunneling times, should pro- 
ceed on similar lines. The precision of the measurement 
scheme is estimated in Sec. [V] 



II. THE DETECTOR MODEL AND ITS 
ARRIVAL-TIME DISTRIBUTION 



We briefly review the detector model introduced in 
Refs. [ill, QjJ [2(| and the arrival-time distribution ob- 
tained from this model in Ref . [l7| ■ The detector consists 
of a three-dimensional array of spins with ferromagnetic 
interaction in a mctastable state. The spins are weakly 
coupled to the environment, modeled as a bath of bosons. 
The effect of the particle to be observed is to strongly 
enhance the spin-bath coupling when the particle's wave 
function overlaps with that of a spin. Thus when the 
particle is close to a spin, this spin flips with strongly 
increased probability due to the enhanced spin-bath cou- 
pling. By means of the ferromagnetic interaction, this in 
turn triggers the subsequent spontaneous flipping of the 
neighboring spins, and finally, by a domino effect, of all 
spins, even in the absence of the particle. In this way, 
the single spin flip induced by the presence of the par- 
ticle is amplified to a macroscopic event, and either the 
change in the detector state or in the bath state can be 
measured. 

The Hamiltonian for this model is 



H — ffpart + #det + -^bath + -f^spon + H coup , (1) 

where we use the following definitions: 

A2 



part 



2m 



(2) 



where a,£ is the annihilation operator for a boson with 
wave vector t 1 is the free Hamiltonian of the environment, 
modeled as a bath of bosons, and 



H spon = ^U j) e if ^ +h . c , 



(6) 



with 



^ id,. (ii, &® = (&v>y = 1^ jiii, (7) 

and the coupling constants 7^ and the phases fl de- 
pending on the particular realization of the detector and 
the bath, describes the permanent spin-bath coupling. 
The spin-bath coupling is strongly enhanced in the par- 
ticle's presence due to 

#c ( £p = e x u) (*) e h + *A , (8) 



with 



U) 



U) 



(9) 



thus, the enhanced coupling of the j th spin to the bath is 
proportional to a sensitivity function x ( x ) which van- 
ishes outside the region where the j th spin is located. 
An example would be the characteristic function which 
is 1 on and zero outside. It will be assumed in the 
sequel that the Markov property [see Eq. (fl6|) ] holds in 
an appropriate continuum limit. Initially, the system is 
prepared in the state 



|*o) = |0)|Ii...|7d)|^o), 



(10) 



where |0) is the ground state of the bath (no bosons 
present), and |^o) denotes the spatial wave function 



of the particle. For the hw^ only slightly above the 
energetic threshold set by the ferromagnetic spin-spin 

( 7 ) 

coupling, and 7^ sufficiently small, the probability of 
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a s pon taneous spin flip (false positive) is very small 
fl8l.ll9l.l20l]. But when the particle is close to the j th spin, 
the excited state ||) • decays much more quickly, due to 

(i) 

the enhanced coupling, g¥ of the spin to the bath. 
Then, the ferromagnetic force experienced by its neigh- 
bors is strongly reduced, and these spins can therefore 
flip rather quickly, even in the absence of the particle, by 
means of the 7^ . The first spin flip will be amplified to 
a macroscopic event by the previously mentioned domino 
effect [lEISISl- 

In Ref. [171 ] we investigated, by means of the quan- 
tum jump approach [2l[ , the application of this detector 
model to arrival-time measurements. The bath modes 
were indexed by wave vectors 



2tt 

Lbath 




m e N , 



(ii) 



and the coupling constants were assumed to be of the 
form 



T? : 

^bath 



-^batli 



gjp = (r«(a,,,e,) + 0(L^ h ))y 
with uji = c(ujg)l, = £/£, and 

r«(^, e ,)| 2 » |rg) on (^,e,) 

We introduce "modified resonance frequencies" 

D 



,~,U) — , ,U) 



E 

fc=i 



(12) 



(13) 



(14) 



which incorporate the ferromagnetic spin-spin coupling. 
The correlation function (t) is defined by 



I 



(J) 



(15) 



(i) (i) (i) 
and similarly for kx' k¥' and ksL It is assumed that 

,37,' 79 ' 77 

in the continuum limit, £bath — *■ oo, the correlation func- 
tions satisfy the Markov property, i.e., 



«SW«0 if r>r c , 



(16) 



for some small correlation time r c , and likewise for the 
other correlation functions. This is the case, e.g., if 
r(u)£,ei) is independent of u>i, as in quantum optics. 

Let IV'cond) denote the time evolution of the spatial 
wave function under the condition that no boson is de- 
tected at least until time t ( "conditional time evolution" ). 
Then 



Po(t) = (V>condkcond) 



(17) 



is the probability that no detection occurred until t, 
which is decreasing in time, and 



(18) 



is the probability density for the first detection to occur 
at time t. In Ref. [T3| it was shown that |V-' CO iid) obeys a 
"conditional Schrodinger equation" of the form 



cond} 



H con d l^cond) 



(19) 



with a Hamiltonian containing a complex potential, 

~ 2 n 



He, 



id = ■z h - [<Wt(x) - iA(£)] , 

2m 2 



(20) 



where <5 s hift(x) and A(x) are defined as follows. A(x) is 
a position-dependent detector decay rate, 

poo 

A(x) = 2Rc£ / dr{4)(T)x«(x) 2 + ^(r)} 
„• Jo 



~U) 



chip 



LOr. C \L) n 



(4^,e)| 2 X «(x) 



r0') (V, (j) e) 

1 spon 1 w i e I 



(21) 



where the df2 e integral is taken over the unit sphere and 
where the contributions from and «^ have been ne- 
glected, due to (fTB"]) . The terms have the familiar form of 
the Einstein coefficients in quantum optics, where there 
would also be a sum over polarizations. The real part of 
the potential, 5 s hift (x), is given by 



4hift (x) = 2 



Im J2 / dr{ 

■ Jo 



4>)x (j) M 2 



(22) 

in quantum optics, this would correspond to a line-shift. 
Since the k 77 term leads to a constant it just gives an 
overall phase factor and can be omitted. 

Note that the Hamiltonian on the r.h.s. of Eq. ([20]) is 
not norm conserving due to the imaginary contribution 
-ihA(x)/2, in accordance with Eq. (J7]). From Eqs. ([20]) 
and l|18p one easily finds 



H, 



cond 



-H. 



t 

cond 



0, 



cond 



d 3 xA(xi) |<x|^* nd>| 



(23) 



which is an average of the position-dependent decay rate 
A(x) of the detector, weighted with the probability den- 
sity for the particle to be at position x and as yet unde- 
tected. 
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It was shown in Ref. [17] , using an example with a 
one-dimensional, simplified detector model, that wi(t) 
essentially agrees with the probability density obtained 
from the discrete model with 40 bath modes by means of 
standard unitary quantum mechanics, up to the time of 
recurrences due to the discrete nature of the bath. 



III. THE RESET OPERATION AFTER THE 
FIRST BOSON DETECTION 

A. The particle reset state 

For the intended application of the detector model we 
will need more than the time evolution up to the first 
detection. In particular, we require knowledge of the 
"reset state" [22|], the particle state immediately after 
the first boson detection. 

Let the complete system (bath, detector, particle) be 
described at a particular time by a density matrix of the 
form 

£? = |o) |Ti • - - To) £> P (Ti • ■ ■ Tz?l (o| , (24) 

(with g p the particle density matrix), i.e., no boson and 
all spins up. If a boson is found in a broadband boson 
measurement a time At later, the density matrix for the 



corresponding subcnscmble is obtained by sandwiching 
the above expression with 

Pi =J2 IW^I, (25) 
i 

by the von Neumann-Liiders projection rule (23l.[24|. and 
the trace gives the probability. The subsystem consisting 
solely of the particle is then described after the detection 
of a boson by a partial trace, 

tr d et tr bath Pj [/(At, 0)q tft (At, 0) Pi = K Qp At , (26) 

which defines the operation 1Z. To calculate this we go 
to the interaction picture w.r.t. 

Hq = H — (ifcoup + H spon ) 

= -ffpart + Hdct + -ffbath 

Hi = H coup + H spon (27) 

and use standard second-order perturbation theory for 
Ui(t, 0) = exp{|i?oO U(t, 0). The zeroth order does not 
contribute to Eq. (f2"6")) . and neither does the first order 
since Pi acts once on |0). In second order only the term 
with Hi on the left and right survives, and in second 
order one obtains after some calculation 



KgpAt = tr dot tr bath P ie -*" 

At 



H At 



-■0 J dtie^ Hotl Hie~i Hotl g ~ J dt 2 ei H ° t2 Hie~^ H ° t2 e^ H ° At Pi 


At ti 

2Ree-i H -' A( ^ / dt x J dt 2 e <^-^) 



,-|ff part At 



(28) 



where x(t) = x + pt/m is the free time development 
of x in the Heisenberg picture. To obtain the second 
equality, the rectangular integration area ti, t 2 € [0, At] 
has been split into two triangles, t\ <G [0, At], t 2 S [0,ti] 
and t 2 g [0,At], ti e [0,t 2 ]. The phases in Hi 
have canceled since we are taking the trace over detector 
states. 



B. An example 

As an example we consider a simplified model with 
only one spatial dimension and only one spin. This sim- 
plification is reasonable if the radius of the region Qj is 
smaller than the distance between spins. (Our assump- 
tion of locality of the interaction though is a bit stronger 
than this however, since below, for calculational conve- 
nience, we will extend the region Qj to a half- line, i.e., 
X(^) - * ®( x )-) The vectors x, p and t are replaced by 
x, p and £, and there is of course no summation over the 
spin index j; we will temporarily drop the superscript . 
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The detector Hamiltonian iJdct [Eq. ([3])] simplifies to 



(29) 



and the modified resonance frequency luq [Eq. (|14[) ] is 
replaced by the resonance frequency of the single spin, 
luq. Also, we will neglect H spon in view of the assumption 

I li\ 2 ^ lff^| 2 - The time development in this model for a 
wave packet incident from the left with detector and bath 
initially in state || 0) has been investigated, among other 
questions, in Ref. (l71 |. 

We assume a maximal boson frequency lo m and N bo- 
son modes, 



cue = uj M n/N, n = 1, 
^ = -iGy/uJt/N. 



(30) 



We further take x(x) = O(x) where denotes Hcavi- 
side's step function. As the particle we consider a cesium 
atom, and the initial state at t = is assumed to be a 
pure state, g p = (ip\, with 



(x\tl>) 



1 



AxV2w 



-x /4(Ai) i(mv /h)x 



(31) 



With these simplifications the first line of Eq. 
yields 



At 



(x \TZg p At\ x) 



* o 

x ^v|0^(*2)]e^ rtAt/fi 



E 



fit 



At 

I dte i{kH - Ua)t (x e - iH ^ {At - t)/h Q {x) e~ iH ^ t/h tpj 



The contribution (x\ - ■ ■ \ ip) has an intuitive explanation: 
The initial state evolves freely until t and is then pro- 
jected onto the detector region. In an intuitive picture, 
the time t may be viewed as time of occurrence of a boson. 
Since a boson can only be created when the spin couples 
to the bath, the particle has to be inside the detector 
at t, hence the projection. Then the state continues to 
evolve freely until At. The integration is understood as 
sum over all possible "paths" satisfying the above pic- 
ture, i.e., as sum over all times t. 

The second line of Eq. (|3"2"]) can be evaluated, e.g., by 
inserting 1 = dk \ k) (k\ with momentum eigenfunc- 

tions (x\k) = i/l/27r e lkx on the left of 0(i), and 

/oo 
dx'dk' \x') (x'\k') (k'\ 
-OO 

T 

2^ 



dx'dk'e ik ' x ' \x') (k'\ (33) 



on its right, and noting that 



(fc 19(5)1*') 



—ikx' 



e(x'). 



(34) 



(32) 



r 



A numerical illustration of (x \TZg p \ x) [note the division 
by At as compared to Eq. (J2HJ) or Eq. 1(52^1] for N = 15 
bosons modes is given in Fig. [5] (dots). 



C. The continuum limit 

We return to the general expression in Eq. ([28]) . For 
simplicity we will assume in the following 



X (j) (x)=x(x), j = !,••• ,D. 



(35) 



In view of Eq. © and the remarks after Eq. (fTU|) we 
neglect the 7^ terms. We introduce the collective cor- 



relation function 



U) 



(36) 



Since we assume that the coupling constants arc such 
that the Markov property holds in the continuum limit, 
i.e., 



k(t) ~ if t > t c 



(37) 
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for some small correlation time r c , in the double integral 
of Eq. (|2"B")) only times with t-y — t 2 < t c contribute, and 
if t c is small enough one can write 



X [x(ti)] « X [x(ta)]. 



(38) 



Then, with a change of variable, r = t\ — 12, the integra- 
tion over r can be extended to oo if t c <C At. Specializing 
the definition Eq. (f2l"j) to the case at hand, we put 



FIG. 2: Dots: reset state (x\1Zq p \x) for a pure state g p — 
\ip) with the Gaussian wave packet of Eq. (|31[1 . Aa; = 
50 nm and w = 1.79 m/s; w = 2.38 x 10 12 s~\ uj m = 
4.6wo, G = 2.782 x 10 3 s" 1/2 , N = 15, and At = 100 uj^ 1 = 
4.185 x 10~ n s. Solid line: reset state A \(x \Q(x)\ %j})\ 2 from 
Eq. (|42|l for the corresponding continuum limit. Up to small 
deviations around x — 0, the reset states from the discrete 
and the continuum model are in very good agreement. 



A = 2Re / dr «(r) 
/o 



(39) 



and note that A is given by the right side of Eq. (j2"Tj) 
without the terms I^pon and x^'(x). One then obtains 



A/ 



R Pp A( = A y , ^ie"* ffpart(A *"* l) x(x)e-^ ff p arttl g p e*^ art * 1 x(x)e^ ffp '' rt(A *"* l) 



A x (x) epX (x) At + 0(At 2 ). 



(40) 



Hence, to first order in At, one finally has 
KQp = ^x(x)e p x(x) 



(41) 



for the state immediately after a detection. 

1Z is called the reset operation. If g p is a pure state, 
then the reset state is also a pure state. In particular, if 

£p = IV'cond) (V'condli tnCn the reSet StatC is S ivCn h V the 

wave function 



preset) = ^4 V2 X(X) |^ „d) 



(42) 



We note that this result is physically very reasonable 
since it means that right after a detection of the particle 
by a detector located in a specific region the particle is 
localized in that region. We also note that [see Eq. (|2"3")) ] 



lib* 

csct V-' reset 



)=w 1 (t), 



(43) 



which is the probability density for the first detection. 

We apply Eq. f4"2")l to the continuum limit of the ex- 
ample of Subsection IIII Bl In that one-dimensional and 



single-spin case one has with Eqs. (j30|) and (|36|) 

N 



«(r) 



| G |2 -iuJ M r (n/N)+iu T 

11 1 L^l AT M 



(44) 



^m|G| : 



J d££e- iu M T t +iuJoT for 



Thus one obtains in the continuum limit 
|G| 2 (l + iw M r)e 



k(t) 



-!(u m -uo)t _ gil^oT" 



A = 2tt|G| 2 — forcj M > w 



<5 s hift = 2 |G|" 



f^ln 


W 


-0 









(45) 



and r c is of the order of ui 1 . The resulting spatial proba- 
bility density A \ (x \Q(x)\ -0) | 2 is plotted in Fig. [2] for the 
same initial state and same parameters as in the discrete 
case. The plots are in very good agreement. 

We note that in the present model with indepen- 
dent of x there is no explicit recoil on the particle from 
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the created boson. This is in line with the original idea of 
a minimally invasive measurement. The absence of an ex- 
plicit recoil distinguishes the present detector model from 
other models which are based on the direct interaction 
with the particle's internal degrees of freedom. An exam- 
ple for such a model would be the detection by means of 
laser induced fluorescence. In that case, the reset state 
after the detection of the first fluorescence photon cx- 
plicitely incorporates a recoil due to the momentum of 
the emitted photon [25[. It appears reasonable that in 
the present model no such recoil on the particle occurs: 
After all, the boson is emitted not by the particle but 
by the spin lattice. Hence, the recoil should be experi- 
enced by this lattice, rather than by the particle, similar 
to what occurs in the Mossbauer effect. Of course, the 
projection of the wave packet onto the detector region by 
means of the reset operation also changes the momentum 
distribution of the wave packet. 



D. Subsequent time development 

After detection of a boson, the further interaction of 
the particle with the detector depends on the particular 
choice of parameters of the detector model. The internal 
dynamics of the detector after the first spin flip has been 
investigated in detail in Refs. [H, HH, H(| • Based on these 
results, several choices are possible such that the ampli- 
fication of the first spin flip will not significantly change 
the spatial wave function after the first spin flip, Preset)- 
This means that effectively only one reset operation, as- 
sociated with the very first spin flip, has to be performed 
on the spatial wave function. Such "minimally invasive" 
detector models will be of interest if one is interested in 
actual quantum mechanical limitations of a passage-time 
measurement. 

As a simple example, consider a ring of identical spins 
with nearest-neighbor interaction, 

utf' = ljq and uj k > = uoj 6 j+ i, k (46) 

with <5j+i,fc Kronccker's S, x^)(x) = x( x ) as before, and 
j = D + 1 identified with j = 1. The rate for the neigh- 
bors of the first-flipped spin to flip into their ground state 
is denoted by A\ (x) and is given by the right side of Eq. 
(|2"Tj) . but with ujq replaced by luq since the ferromagnetic 
forces on these neighboring spins cancel. Choosing pa- 
rameters as, e.g., loq 3> ujq and T, r spon independent of 
lu£ as well as c{u>t) = Co, one has A\ 3> A. By a kind 
of domino effect the whole ring flips into the ground 
state; the mean time needed for this given by D/2Ai 
El El III- H this time is very short, as one can achieve 
by making loq large (while loq remains small to prevent 
spontaneous spin flips before the first particle-induced 
spin flip), the reset operations associated with these sub- 
sequent spin flips will not significantly change the par- 
ticle's state since the wave function has been projected 
onto the detector by the first reset operation already. 



Another possibility to prevent the spatial wave func- 
tion from being changed by the amplification process 
would be to couple only one spin to the particle, by 
choosing = geSjj in Eq. ©. (In fact, this is the de 
fac to setup of the detector actually investigated in Refs. 
HE EE H0| : Effectively only one spin couples to the par- 
ticle, and subsequently the other spins flip spontaneously, 
i.e., without particle-enhanced spin-bath coupling.) 

IV. APPLICATION TO PASSAGE TIMES 

A. General setup 

We now consider two detectors separated by some dis- 
tance. As indicated in the preceding section, we assume 
that the amplification of the first spin flip to a macro- 
scopic event is very fast and does not change the spatial 
wave function. Thus, we take the probability density 
wi (t) [Eq. ([23]) ] for the first spin flip to be the "measured 
arrival-time distribution" and Preset) P^l- (H2])] as the 
particle's state after the detection. Then, the joint prob- 
ability density for the first detector detecting the particle 
at T and the second one detecting it at T + r is given by 

G(T,T+r) = i»W (T; |^» (r; ) , (47) 

V || I Preset) 11/ 

where the superscripts indicate the detector under con- 
sideration and the second argument is the initial state 
for which the respective probability density is calculated. 
Since wx is bilinear in the wave function [see Eq. (J23J)], 
this simplifies to 

G(T,T + r)=^ 2) (r;|^ T CS ct)) (48) 

by Eq. ([43]) . The desired measured passage-time distri- 
bution is then given by integration over the entry time 
T: 

G(r) = J dTw™ (r;|i&Lrt»- (49) 

B. Numerical example 

In order to investigate basic features of a quantum 
passage-time distribution obtained from the detector 
model in the continuum limit and in the above mea- 
surement scheme, we consider as a simple example a ce- 
sium atom in one dimension. The initial wave packet 
is prepared in the remote past far away from the detec- 
tor such that the free wave packet (with no detectors 
present) would be described at t = by a Gaussian min- 
imal uncertainty packet around x = with Ax = 1/im 
and average velocity vo = 0.717cm/s. Each of the two 
identical detectors is described in the continuum limit by 
an absorbing potential — iV = —ihA/2 extending from 
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to 20/im or lOO^m to 120/im, resp., where we consider 
three different examples A = 1.4337 x 10 3 s"\ 2.3895 x 
10 3 s^ 1 , 2.3895 x 10 4 s" 1 . These parameters are chosen 
in such a way that transmission and reflection without 
detection, which are typical for imaginary potentials [26| 
and have been extensively studied in the framework of 
quantum arrival times [27], HH, [291 . [30j . play no signifi- 
cant role. Consequently, all distributions shown in the 
following arc normalized to 1 to a good approximation. 
The passage-time distributions, calculated as described 
in Section IIV Al are shown in Fig. [3] 
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FIG. 3: Passage-time distributions calculated from the detec- 
tor model in the continuum limit for three different values of 
the rate A for the first spin flip in the presence of the particle, 
all other parameters kept fixed; for comparison, the dotted 
line shows the passage-time distribution for an ensemble of 
classical particles which have the same momentum distribu- 
tion as the initial wave packet. 



It is seen that small and large values of A give rise to 
rather broad distributions, while the intermediate value 
yields a narrower one. The reason for the broad distri- 
bution arising for small A can be understood by looking 
at the arrival-time distribution measured by the first de- 
tector (see Fig. 01): It is already this distribution which is 
rather broad for small A. Physically, small A means that 
the detector is responding only slowly to the presence 
of the particle; the undetected amplitude |Vcond) decays 
only slowly, yielding a broad arrival-time distribution 
at each of the two detectors and consequently a broad 
passage-time distribution. In other words, the poor qual- 
ity of the passage-time measurement arises from the poor 
quality of the individual arrival-time measurements. The 
measurement of the arrival time, however, can be im- 
proved by making A large, that is, by making the de- 
tector responding faster to the presence of the particle. 
In Fig. Hit is seen, e.g., that for A = 2.3895 x 10 4 s" 1 
one comes much closer to Kijowski's arrival-time distri- 
bution; this in turn is known to have minimum standard 
deviation among all those distributions fulfilling certain 
axioms transferred to quantum mechanics from classical 
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FIG. 4: Arrival-time distributions as obtained from the first 
detector. The detector with small A gives a broad distribu- 
tion indicating poor quality of this measurement. Enlarging 
A carefully (to avoid reflection), one may approach Kijowski's 
axiomatically derived arrival-time distribution at x = (dot- 
ted line). 



arrival-time distributions [3l| . and thus provides a sort 
of "ideal distribution." 

While the broad passage-time distribution for small 
A can be understood simply as due to the poor qual- 
ity of the individual measurements, the broad distribu- 
tion for large A is more interesting. It can be under- 
stood looking at the reset state immediately after the 
detection by the first detector. Of course, this reset 
state depends on the instant of time when the detection 
took place; as an example, we consider detection times 
Ti = 0.041 ms if A = 2.3895 x 10 4 s" 1 and T 2 = 0.167 ms 
if A = 1.4337 x 10 3 s _1 , which are close to the maximum 
of the respective probability distribution w\(t) (see Fig. 
U). The reset states and the free packet at the respective 
times are shown in Fig. [5] It is seen that the fast de- 
tection has a strong impact on the wave function. Large 
A means that IV'cond)' m particular that part of IV'cond) 
which overlaps with the detector, decays very fast. But 
since the reset state (|42l) immediately after the detection 
is essentially the projection of IV'cond) on *° * ne detector 
region, the fast decay of this overlap yields a reset state 
which is very narrow in position space, located at the 
very beginning of the detector. Thus, by the Heisen- 
berg uncertainty relation, it is very broad in momentum 
space as can be seen in Fig. [5] It is intuitively clear that 
such a broad momentum distribution immediately after 
the measurement of the entry time by the first detector 
yields a broad passage-time distribution. So the broad 
passage-time distribution in case of large A is due to the 
strong distortion of the wave function by the first mea- 
surement. Since the broad momentum distribution of the 
typical reset states in this case is enforced by the Hciscn- 
berg uncertainty relation, this is a pure quantum effect. 
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FIG. 5: The normalized reset state in position space im- 
mediately after a detection at (a) T2 = 0.167 ms (with 
A = 1.4337 x 10 3 s" 1 ) and at (b) Ti = 0.041 ms (with 
A — 2.3895 x 10 4 s _1 ), compared with the free wave packet at 
the respective instance of time (dotted line). The fast detec- 
tion in case (b) has a strong impact on the wave function. 
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FIG. 6: Momentum distributions of reset states immediately 
after a detection at T2 = 0.167 ms (solid line, with A = 
1.4337 x 10 3 s" 1 ) and at Ti = 0.041 ms (dash-dotted line, 
with A = 2.3895 x 10 4 s _1 ), compared to the momentum 
distribution of the initial wave packet (dotted line). The fast 
detection in the latter case leads to a strong broadening of 
the momentum distribution. 



V. WIDTH OF THE PASSAGE-TIME DENSITY 

In Rcfs. [13, EH a measurement scheme for passage 
times was investigated, based on the (semi-) continuous 
coupling of a particle to a clock. It was argued in these 
references that the precision of the measurement behaves 
like -E -1 where E denotes the kinetic energy of the par- 
ticle. One now may wonder whether or not the E~ x 
bchavior is a fundamental quantum limit for measuring 



passage times. We argue that this is not the case since 
the present measurement scheme by means of two spa- 
tially separated detectors yields, for optimal parameter 
choices, passage-time densities with widths behaving like 
-E~ 3 / 4 . Thus, for low energies, we have an example which 
breaks the E~ x limitation of the clock model. 



A. Estimating the precision 

In this subsection we give an estimate for the width of 
the passage-time density obtained from the present mea- 
surement scheme. We assume that the detectors can be 
described by the continuum limit and that transmission 
and reflection without detection are negligible. This as- 
sumption is justified in the examples of the preceding sec- 
tion, which employed rectangular sensitivity functions. It 
can also be justified in general if one drops the restriction 
to rectangular sensitivity functions xi x ) 

Considering particles with mean velocity vq, the detec- 
tor is assumed to be constructed in such a way that the 
first detection occurs with high probability in a spatial 
interval of length L. The length L is related to A of Eq. 
(f59")) . an average detection rate, L being of the order of 
vq I A. The length L imposes an upper limit on the width 
of the reset state in position space, 

Ax rcsct < L w v /A. (50) 

By the Heisenberg uncertainty relation, this immediately 
yields a lower bound for the width Ap rcsot of the reset 
state in momentum space, 

Ap rcsot > h/2Ax ICBCt > h/2L. (51) 

Note that this is only a very rough estimate, without 
taking into account details of the actual wave packet. 
If the incident wave packet is very narrow in position 
space, then Ax rese t also may be much smaller than L, 
and consequently Ap reset may be much larger than H/2L. 
Also, the reset state may be far from being a Gaussian, 
and then already the first inequality in Eq. (|5i"j) may 
underestimate the width of Ap rcsot seriously. 

Let At denote the width of the measured passage- 
time distribution. There are several contributions to this 
width. First, a particle with velocity vq takes the time 
r = d/vo to travel the distance d between the two detec- 
tors, and therefore the width of the reset state in position 
space contributes 

ATreset.K = Aa; reS et /^O (52) 

to At. Second, the width in momentum space con- 
tributes according to 

_ Ap rcsct d h d 

"Preset ,p — _ (~. \ n ■ {^^) 

mv Vq 2 Aa; reset UIVq 

Third, the width of the delay, Ardeiay, for the first detec- 
tion due to the spin-boson interaction is roughly 

Ar dolay = 1/A w L/v . (54) 
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An estimate for the width of the passage-time distribu- 
tion is given by the sum of these contributions, where 
Aideiay has to be counted twice since it arises in both 
detectors: 

At — ^ATdelay -|- At ycsc ^_ x -\- AT Yese i^p 

> 2/ A + Ax Teset /vo + M/2mv 2 ,Ax rosQt . (55) 

From this estimate it is again seen that both small as well 
as large values of A, i.e., both slow as well as fast detec- 
tors, lead to rather broad passage-time distributions (due 
to A-r dG i ay ~ 1/A and Ar rese t,p ~ 1/Ax resct ~ 1/L ~ A, 
respectively). 

B. Optimal parameters 

Having established the general estimate for At in Eq. 
(|55p . we now turn to the task of finding optimal parame- 
ters, minimizing At. We are interested in measuring the 
passage time through a spatial interval of length d, the 
distance between the starting points of the two detectors, 
which we regard as fixed. First, we consider given detec- 
tors, i.e., a given detection rate, A. Again, particles with 
mean velocity vq would be detected within an interval 
with length L given approximately by L w vq/A. This 
means that the velocities must not be too large since in 
order to avoid undetected transmission L must not ex- 
ceed the actual length of the detector (and in particular 
L must not exceed d). Quantum effects, however, are 
expected to play a role for slow particles while fast parti- 
cles may be treated classically, hence this is not a serious 
drawback. 

We assume that the reset state is a Gaussian wave 
packet, or at least close to a Gaussian, so that in the 
second line of Eq. (|55p the approximate equality holds, 

At r* 2/ A + Ax rosot /v + M/2mvlAx rcsct . (56) 

We note that 2/ A is a purely detector-related quantity, 
while the remaining contribution to At is determined by 
the shape of the reset state only. Wc will first optimize 
this latter quantity. For given particles with given mean 
velocity, i.e., m and vq fixed, this is minimal for 

Preset = V M / 2mv 0- (57) 

Substituting this into Eq. (|56p and writing E = mv^/2 
for the kinetic energy of the incident particles yields 

AT opt; rosot « 2j A + y/h^/mj^ £~ 3/4 . (58) 

The subscript "opt; reset" indicates that only the reset 
state was optimized while the detection rate A was con- 
sidered as a given quantity. 

Aiming at optimizing also the detection rate A for min- 
imal At, one would like to choose A as large as possible 
in order to reduce the 2/ A contribution to At. How- 
ever, one has to take into account that the width of the 



reset state is bounded by Eq. (|50|) . Thus, given d, m, 
and ^o, the decay rate must be at most of the order of 
vq I Aa^gget with Ai^ot as in Eq. ([57]) . In fact, we may 
choose the incoming state such that it forms a Gaussian 
minimal uncertainty packet at the starting point of the 
first detector with width in position space 

Ax = Ax°ps Ct = \Jhd/2mvQ, (59) 

and choose further 

A « v a /2Ax = y/m/2M v 3 / 2 . (60) 

We note that, by this parameter choice, yet another re- 
quirement of a good measurement scheme is fulfilled: The 
detection by the first detector will not change the wave 
function too strongly. The reset operation after this de- 
tection is essentially a projection onto the detector re- 
gion, and the detection is slow enough that at typical de- 
tection times most of the wave packet overlaps with the 
detector, hence the projection does not change the wave 
function too much. Thus, the reset state will be close 
to a Gaussian wave packet with width Ax = Ax°^ set . 
Substituting Eq. (SO) into Eq. (55) finally yields 

AT opt » E- 3 /\ (61) 

We stress that, independent of the detection rate A, 
the optimal energy dependence of At is limited to E~ 3 / 4 
already by the dependence of At on the width of the reset 
state in position space, AT rese t x ~ Ax reset [Eq. (|52p], and 
on its width in momentum space, Ar rese t iP ~ Ap reset ~ 
l/Ax rose t [see Eq. (53)]. 

For the example of a cesium atom with vq = 0.717 cm/s 
and a distance of d = 100 /jm between the detectors (with 
rectangular sensitivity function) of the preceding section, 
optimal values would be according to Eqs. ([57]) and ([60]) 

Aa; pt = 1.83 fxm, and A opt = 1.959 x 10 3 s" 1 . (62) 

Considering the wave packet with Ax = 1 /urn actu- 
ally investigated in the examples, the optimal decay rate 
according to Eq. ([6"0]) would be A opt (Ax = 1 /Ltm) = 
3.585 x 10 3 yus -1 . This is consistent with the observation 
that the example with A closest to A opt (Ax = 1 /Ltm) 
yields the narrowest distribution. 

Summary 

We have investigated the continuum limit of a fully 
quantum mechanical spin-boson model for the detection 
of a moving particle and its application to passage-time 
measurements. The continuum limit has been derived 
under the condition that the spin-boson interaction sat- 
isfies the Markov property an assumption that was ex- 
plicitly verified in Ref. [17j. Analytical expressions for 
the state immediately after the first detection have been 
obtained, and for an example with a simplified detec- 
tor model and 15 boson modes it was shown numerically 
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that the continuum limit is in very good agreement with 
the discrete model. Further, analytical expressions for 
the passage-time distribution have been obtained, and 
numerical examples for passage-time measurements have 
been discussed. Detectors with a very slow response yield 
broad distributions, due to the poor quality of the indi- 
vidual measurements, and so do very fast detectors, due 
to the strong distortion of the particle's wave function 
by the measurement. Intermediate detectors yield nar- 
rower distributions. The optimal precision of the present 
measurement scheme has been estimated to behave like 



E~ 3 / 4 , where E is the kinetic energy of the incident parti- 
cle. For slow particles this is better than, and in contrast 
to, a scheme based on quantum clocks which yields an 
-E -1 behavior. 
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